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£f} ' Abstract. An elliptic analogue of the q deformed Knizhnik-Zamolodchikov 

equations is introduced. A solution is given in the form of a Jackson-type 
integral of Bethe vectors of the XYZ-type spin chains. 

> , 

CN , Introduction 

o 

■ In this paper we introduce a holonomic system of difference equations associated 

to elliptic R matrices and give its solution in the form of Jackson-type integral, 
following Reshetikhin's idea Q for the trigonometric R matrices. 

Reshetikhin constructed a solution to the q-deformed Knizhnik-Zamolodchikov 



equations FR] by a Jackson- type integration of Bethe vectors of the XXZ-type spin 
chain models. Matsuo [Ma] also found the same kind of formulae from a different 



viewpoint. 

q-( On the other hand, the Bethe Ansatz method for the spin chain models associated 

to the elliptic R matrices has been studied since Baxter Q. Hence a natural 
question is how to find an elliptic version of Reshetikhin's approach to the q-KZ 
equation. It turns out that the argument in B can be carried out for the elliptic 



R matrices as well, except for one point. In contrast to the trigonometric case, 
an elliptic spin chain model does not have a unique vacuum vector in its local 
state space but a series of "pseudo- vacua" which depend non-trivially on a spectral 
parameter. This dependence breaks down naive analogy. 

We overcome this difficulty by introducing a "space of Bethe vectors" and a 
boundary operator which shifts a spectral parameter. 

In the vertex picture the linear space of Bethe vectors depends on spectral pa- 
rameters. Therefore we have to use the IRF picture in order to interpret the system 
as a holonomic matrix difference system in the sense of Aomoto Q. In this con- 
text, our system is described in terms of representation of Felder's elliptic quantum 
groups fj. 

This paper is organized as follows. In the first section we recall several facts 
related to the elliptic R matrices and introduce a space of Bethe vectors and a 
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boundary operator. We define a system of difference equations in the next section 
and show its holonomicity. The third section is devoted to construction of a solution 
of this system by a Jackson-type integral of Bethe vectors. 

Acknowledgements. The author expresses his gratitude to Nicolai Reshetikhin 
for comments and discussions, and to the Department of Mathematics of the Uni- 
versity of California at Berkeley for hospitality. He is supported by Postdoctoral 
Fellowship for Research abroad of Japan Society for the Promotion of Science. 

1. Space of Bethe vectors 

In this section we define a space of Bethe vectors of an elliptic spin chain and 
linear operators acting on this space. The state space Tt of a finite XYZ-type spin 
chain is defined to be a tensor product of the local state spaces: 

(1.1) H := V h ®V h ® ■■■®V lN , 

where V lj (j — 1, . . . , N) are the spin lj representation spaces of the Sklyanin alge- 
bra. (See Appendix |a| for a review of the Sklyanin algebra and its representations.) 
We will introduce a direct sum of subspaces of TL depending on spectral parameters 
which we will call a space of Bethe vectors. On this subspace act not only the R 
matrices but also a boundary operator which shifts the spectral parameters. 

We assume that the parameter r\ which determines the structure constants of 
the Sklyanin algebra is a rational number, r\ = r' jr. This assumption is necessary 
to consider analytic solutions of the system later. 

Space of Bethe vectors. Baxter Q introduced vectors which intertwine vertex- 
type Boltzmann weight of the eight- vertex model and IRF-typc Boltzmann weight. 
They are generalized to higher spin cases, (cf. DJKMOj, fTlf)) 



Definition 1.1. An (outgoing) intertwining vector, <f>\ \i (it), is a vector in the spin 

I representation space V 1 = Oqo + of the Sklyanin algebra, defined by: 

(1.2) 



+m / A — \ / A — \ 

II *io [v + + (2j - 1 l)v) Oio [y-^-{2j-l- 1)„J 



3=1 

x II (v + + ^ - 1 - V*) 9w ( y - ^r 1 - ( 2j 

3=1 

where A, A' are parameters satisfying A — A' = Amrj, m G {—I, — I + 1, ■ • • , I}, and 
u is a complex parameter called a spectral parameter. We call 

(1-3) 4°(«) = *SxW«) 

a local pseudo-vacuum. 

It is easy to see that generically {</>a+4«"?, \( u )}m=-i, -i+i,. .. ,i is a basis of V ■ 
Graphically an intertwining vector is denoted as in Figure [l]. 



1 The normalization in previous papers [Tl by the author has been improved. Several com- 
plicated factors are absent now. Vectors <f>^ x' ( u ) m those papers correspond to the following 
4>y a( w ) m * ne present paper. 
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Figure 1. An outgoing intertwining vector <f>\\t(u). 
Definition 1.2. A path vector is a vector in H defined as follows: 

o 

(1.4) \a , ai,... , a N ; zi, . . . , z N ; A) := 

where ao, Oi, . . . , ajy are integers, satisfying the admissibility condition, 

(1.5) a 3 -_i - a 3 G {— 2Zj, — + 2, . . . ,2Z i - 2, 2/j}, 

o o 

Zj's are complex parameters, A„ = A + 2a n rj and A is a fixed parameter which we 
omit hereafter unless necessary. We call a path vector 

Qh<-> l N(zi,... ,z N ) = \a ,ai,... ,a N ;zi,... ,z N ) 

!-6) 

= ^xM) ® w£ (z 2 ) ® • • • <8 

o 

with a = a, a,- = a + 2(Ji + • • • + Aj = A + 2a,?7, a (global) pseudo-vacuum. 

We often denote fi„ '"' ' 1n (zi, ... , z N ) by Sl a (zi, ■ ■ ■ > %n) f° r simplicity. 
Since n = r' jr is a rational number, 

(1.7) 

I ao + r,oi +r, ... ,aiv + r; zi, ... ,zat) = |a ,ai, ... ,ojv;«i, • •■ ,^jv), 

because of the periodicity of theta functions. 

A graphical notation for a path vector |<zo, Oi, . . . , ojv; Zi, . . . ,Zn) is shown in 
Figure f§ Though we should write A„ in this figure instead of a n for the consistency 
with Figure EL we use a n for later convenience. 
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Figure 2. Path vector \clq, Oi, . . . , ajy; zi, . . . , zjv). 




Definition 1.3. The space of Bethe vectors 03^' with spectral parameters 
z = (z%, Z2, ■ ■ ■ , Zn) is a space of functions 

/ : Z/rZ 3 is i ► /(i/) E H 

of the form 

r-l 

fi v ) = X! e 27 " 10 ^ ^2 \ a o, a!,... , a N ; zy, . . . , z N ), 



a=0 



where <ps are complex numbers and sequences of integers a — (oq , a\ , 



satisfy cto = ajy = a and the admissibility condition (1.5). (cf. (1-7).) 

Generically, 
(1.8) 



e 'ao, 



, z N ) 



fflo = «]v = ae{0,...,r-l}, 
(ao, . . . , qn) satisfies (1.5). 



is a basis of 03' 
(1.9) 



Zl,... ,Z]y ' 

dimQ5lV- ,Z . 



Thus 

f = r dim(weight zero space of W lx 



■W 



where W are the spin I (i.e., (21 + l)-dimensional) irreducible representations of 
the Lie algebra sZ(2,C). 

Remark 1.4. When r\ is not a rational number (or, exactly speaking, not a point of 
finite order on the elliptic curve C/Z + rZ), the sum should be taken over all a G Z 
and v is a continuous parameter. The sum might be considered as a formal series. 

R matrices. Elliptic R matrices R = R 1 ' 1 (it) acting on the space V 1 ® V 1 are 
constructed by means of the fusion procedure, (cf. ]Ch| , |DJKMO|, HZj ], [ T2 ] . ) 
We recall the following most important properties and refer details to [T2fl. 

(i) R 1,1 (u) is a linear endomorphism of V 1 ® V meromorphically depending on 
a complex parameter u. (Figure [|) 

(ii) ( Yang-Baxter equation) 

(1.10) R l { 2 (zi - z 2 )R l {l (zi - z 3 )R l 2 ^ (z 2 - z 3 ) = 



R 



i', i" 

23 



(z 2 - z 3 )R l / 3 (zi - z 3 )R l { 2 (zi - z 2 ) 



1,1' I 



as an endomorphism of V 1 ® V 1 ®V l . (Figure 
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Figure 4. Yang-Baxter equation. 
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Figure 5. Unitarity. 



(iii) (Unitarity) 
(1.11) 



R l {\ (u-v)R l 2 { l (v-u) 



Id, 



V l ®V 1 ' ' 

as an endomorphism of V 1 ® V 1 . (Figure ||.) 
(iv) R 1,1 (u) acts on the intertwining vectors as follows: 

(1.12) & l \u-v)<j>®(u)®<j>$Jv) = 



W 






A' 


u — v*J 









where the sum on the right hand side is taken over \j! satisfying fjf — fj, = Amq 
(m G {-/, -l+l, . . ■ , 0) and X—fx' = Am!r) (m G {-I', -I'+l, ... , I'}). Scalar 
factors W are the Boltzmann weights of the IRF-type model. (Figure ^. In 
the figure W is denoted by a crossing of dashed lines.) 

Thanks to (iv), the R matrix defines a map between the spaces of Bethe vectors: 



(1.13) iWfe 



Zj+l 



'■- P jj+lRj,j+l( z j ~ z j+l) : %$zu 



,lj +1 ,lj,... ,1? 
,Zj +1 ,Zj,... 
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Figure 6. R and IRF weight are intertwined. 



where Pj,j+i is a permutation operator of the j-th and the (J + l)-s t component of 
the tensor product V l1 ® ••• ® V N . With respect to the basis (1.8), the R matrix 
is described by the IRF- type Boltzmann weight, and thus by a representation of 
Felder's elliptic quantum group Jf]]. For our purpose, the following special cases of 
(1.12) are important. 

When both 4>\,\'{ v ) an d 4>\'\( u ) are local pseudo- vacua (|Q|), i.e., when A' = 
A + Alrj and /i = A' + Al'r], equation (1.12) is simply 

(1.14) R l > l '(u - v)4\u)^4\(v) = 4l illv (u)®4'Hv). 



When / = 1/2, the R matrix R 1 / 2 ' 1 is expressed as the L operator ( A..4) through 
the identification (A. 9): 



(1.15) 

pi/a,!, x Mgg) i (T ( , ^ f WwiM (0 

* {U) = U« + (2/ + 1),/ + ??)) = ^ (? n (« + (2Z + 1)„) ° ® P {S } ' 



a=0 



In particular, i? 1 / 2,1 / 2 (u) is proportional to Baxter's R matrix ( |A.5| ). Likewise the 
intertwining vectors 4>^x±2„( u ) correspond to the following vectors in C 2 

(116) 6 {1,2) (u-n)-c(- e ° l{iX±u)/2;T/2) 

{1Ab) <P\±2v^ u T1 '- C { 9 0Q ((X±u)/2;t/2) 

where C is a constant: 

M0;t) 2 M0;t)M0;t) 



C = e _7rir / 8 



26 w (0; r/2)0 O i(O; 2r)0 lo ((l + t)B 10 {(1 - r)/4; r) 



Under this identification, the action of R 1 / 2 ' 1 (1.12) can be stated in the following 
form which will be used later. Let us define a matrix of the gauge transformation 

by 



(1.17) 

'-0„i((A-u)/2;t/2) -0„i ((A + «)/2;t/2)\ (l 
9oo((A-u)/2;t/2) 9 qo ((A + «)/2;t/2) )\ ^i(A;r)- 



Af A («) := 
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and a twisted L operator by 
(1.18) 

'ax,x'{u\v) P\,\>{u;v) 



j x,x 



(u;v) 



M x (u)- 1 L{u-v)M x ,(u). 



Then the matrix elements of L x ,x'(u;v) act on an intertwining vector as follows: 
Let A - A' = Arnr) (m £ {-/, -l+l,... , I}). Then, 

(1.19) 

9n(u -v + 2mr])6 11 (X + 2(1 - m)r)) m 



ax,x'(u;v)(j) x l ) x (v) 



hi(u — v + A — 2mri)0ii(2{l + m)r}) ji) 



A'-2?j,A-2?j 



en(A)fli 1 (A-4m» 7 )fli 1 (2»?) 
?h(m - i; - A + 2mrj)9 n (2(-l + m)ry) (i ) 



611(277) 



A'+2jj,A-2jj( W )' 
A'-2)) : A+2jj( W )' 



(1.20) 

/?a,a'( u ; w )^a' ) ,a( w 

(1.21) 
(1.22) 

%(A- 4^)^(2,) 

This formulation was found in the context of the quantum inverse scattering method 
and the algebraic Bethe Ansatz |TF]. See also [Tl|. Note that the column vectors 
of M x (u) are essentially outgoing intertwining vectors. Therefore, defining the 
incoming intertwining vectors 4>\.\>(u) as row vectors of M\(u) . 



hi{u -v - 2mr))0u(\ - 2(1 + m)rj) (;) 

< P\'+2r 1 ,\+2r 1 \ V )- 



(1.23) 



M x (u)- 



/ 0A,A-2 I? (m - TlY 

^x,x+2 V (u - rf) / 

and denoting them as in Figure ^, we can rewrite formulae (1.19)— (1.22) as in Fig- 
ure H (Exactly speaking, the normalization here is different from that in Figure ^, 
which is not essential. In general, incoming intertwining vectors are defined as a 
dual basis of {(j)x+imrt,x{u)} m =-i-i+i, See B.) 



A' 



Figure 7. An incoming intertwining vector </>a' iA i(m)- 

Especially a, 7 and <5 act on a local pseudo vacuum as follows: 

(1.24) ax,x>(u;v)u{ l )(v) = a l (u - v^l^v), 

(1.25) 1 x,x>(u;v)J x l )(v)=0 

(1.26) Sx,x'(u;v)u{ l )(v) = S l (u ~ v)J$ +2tl (v), 
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U — 7] 










A 



a*' ; A 

u — r\ 



Figure 8. An element of a twisted L operator acts on <fi\>,\(v). 



where 
(1.27) 



a\u) 



[u + 2ln) 



l(2t?) 



5\u) 



6 n {2r)) • 



Boundary operator. We introduce an operator Z which shifts a spectral param- 
eter of a Bethe vector. Let k and c be fixed complex parameters. 

Definition 1.5. A boundary operator Z = Z^"' ' l *L (k, c) is a linear map from 
tc > defined by 



(1.28) Z( e 2 ™^>o, • ■ ■ , ajv;Zl) 



e 2 7rl a J v-i^ e c(a J v-a J v- 1 -2/ JV ) ^ 



, aN-V,ZN + K, 21, . 



, 2AT_i> 



(Note that oq — by Definition 1.3. cf. Figure KJ) 



With respect to the basis (1.8 ), Z is a composite of a permutation and a diagonal 
matrix. 

The following property is important. 
Lemma 1.6. The operator Z commutes with an R matrix as follows: 
(1.29) R^-^izN-i z^Z^^^^k^Z^-^c) = 

ry^N — l, il,-" i^N— 2,In ( ,„ J\ ryU ,. . . 1 N _ 2 ,l TV ,In -1 ( ,~ „\ D^ - l^N 



2iH 

*il i 



In — 2,1 n Jn- 

Z\,... ,Zjy-2,ZN ,ZN 



as a map from to fB^^X^HU- 



This is easily proved by comparing the action of the both hand sides on the basis 
( ITsl ) . A graphical notation for this equality is shown in Figure [l(]. 



2. Difference equation 
In this section we introduce a system of difference equations which is an elliptic 



analogue of the equations introduced in |FR|, and show holonomicity of this system 
in the sense of Aomoto 01 . 
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Figure 9. A boundary operator. 



ZN-l 



a 



ai 



1»-2 I O.N-1 i aN 




0- 



ZN + K 



R 



ZN-l + K 



Zl 



ZN 
ZN-l 



a 



«i 



i a/v-i i aN 



zn-i + k 





®Z N + K 
1 



Zl 



R 



z N . 



ZN 



Figure 10. RZZ = ZZR, Lemma 1.6. 



Let us first define a linear operator Aj(z), z — (z\, . . . , Zn), from S^ 1 ' ;; to 



23 



Zi,... ,In 

Zl,... ,Zj-\-K,... ,ZN 



by 



(2.1) Aj (z) := k\ V/ (*j + « - • ■ ■ i?^" 1 + « - zi) 



T. 



z l 1 ,...,l j - 1 ,lj +1 ,...,l N ,l j ( r \phM ( Z .- Z . T ) E 3,j+1 (z--z-^\ 
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or, equivalently, 

(2.2) Aj (z) := R^jlj 1 (zj +k- zj-i) . . . R^ 1 (Zj +k-Zi) 

Z«^ K ,c)R l i^(z j -z N )...R l ifiX 1 (z j -z j+1 ^ 

where 

(2.3) Z»( K) c ) = Pj , ; /' ; 2J , ...P 12 



rt 



dj-1 



a j+ i 



a N 



■© 



Zl Zj-1 Zj + K Zj+l Z N 

Figure 1 1 . Operator Aj (z) 

Then we can define our main object, a system of difference equations: 
(2.4) f(z j )=A j (z)f(z), 



, N where Zi 



for j = 1, 
panding / as in (f.3), 



, Zj + K, . 



z N ) and f(z) e *8%;;;:%. Ex- 



/(£) = ]>> 2 ™]>>(z)|a;z}, 



a=0 



we can regard the system ( [2.4| ) as a system of equations for fg(z). In this IRF 
picture, this system is holonomic in the sense of Aomoto Ja|] due to the following: 

Proposition 2.1. Operators Aj{z) are compatible: 

A (z k )A k {z) = A k (zj)Aj(z). 



This follows from the Yang-Baxter equation (1.10), the unitarity (1.11) and the 
commutativity of R and Z (Lemma 1.6) as in Theorem 5.4 of [ FR] . Symbolically, we 
have only to change the order of crossings of lines in Figure ^ , using the procedures 
Figure ||, Figure || and Figure [l0[ 



3. Jackson-type integral solution 



In this section we give a solution of the system of difference equations (2.4) which 
is expressed as a Jackson-type integral of Bethe vectors. This result is an elliptic 
analogue of that of Jr]. We assume that M := h + ■ ■ ■ + In is an integer. 
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Figure 12. Holonomicity of the system (2.4). 



z N 



' — N 



First let us define the monodromy matrix T(u; zi, . . . , z/v) by 



(3.1) T<~ — ^{u;z N ,... ,zi) 



a^'—'^u;^,... ,21) D 1 "'—^;^,... 

:= L l $ ( u -z N )--- 4 2 (u - z 2 )L\ [u - z x ) 



as an endomorphism of C 2 ®7i = C 2 ® y' 1 ® • ■ 
only on C 2 and : 



>V ij *. ReveL J (u) acts non-trivially 



a=0 



Pi ' =1 



1 O <8> 1 • • • <E> 1 : t/ T)7) (sZ(2)) ~> End(ft). 
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Bethe vectors of the XYZ type spin chains are constructed by means of the twisted 
monodromy matrix defined as follows: 

(3.2) T$;- h (u;z N ,...,z 1 ) = 

_ f Al \,X'" ' h ( u > Z N, ■ ■ ■ > *i) B l ^y-' h {u;z N , . . . ,zx)\ 
^(7jjv,.-,'i^. 2jV) D l ^,-' h (u;z N ,... ,zi)J ' 

:= Mxiu)- 1 ^'- > h (u; z N ,..., zi)M v (u), 



where M\{u) is defined by ( 1.17 ). We often denote, for example, B N '"' ' \ by 

\+2ar),\+2a'T] 

B^a' "' 1 or B a ,a' f° r simplicity, and graphically as in Figure [l3| (cf. Figure ||.) It 



a' + 1 



Zl 



ZN 



a - I 



Figure 13. Operator B l *£" ' h (u; z N , .. . ,zi). 
follows from (1.24), (1.25), (1.26) that 

(3.3) A l ^l(u; z N ,..., Zl )n a (z) = (j[ a l > (u - Zj )\ ^(z), 

(3.4) C^ 2 ^(u;z N ,... , Zl )Q a (z)=0 

(3.5) D l £ 2Ai h a (u; z N ,..., Zl )n a (z) = I J[ (u - z s )\ n a+1 (z). 



By a standard argument in [ TF ] , these operators satisfy the following commutation 
relations: 

B a>a > + i(u;z N ,... ,Zi)B a+ i ta >(v;z N ,... ,z ± ) = 

=B a<a > +1 {v\z N ,... ,zi)B a+l!a >(u;z N ,... ,zi) 
A ata ' + i(u;z N ,... ,Zi)B a+ i ta ,(v;z N ,... ,Z\) = 

=a(u - v)B ata ^ 1 (v;z N , ... , zi)A a+1 ^(u]z N , ... ,Zi) + 
+f3 a >(u - v)B a . a ,^ 1 (u; z N , . . . ,zi)A a+1>al (v;z N ,... ,Zi), 
fla-M'(«;zjv,-- - ,Zi)B a , a >-i{v\z N , ... ,Zi) = 

=a(v - u)B a+ x ta ,(v; z N , . . . , zi)Z? QjQ /_i(u; z N , . . . ,zi) - 
-P a (u - v)B a+lia >(u; z N , . . . , zx)D a ^^x(v; Z N , . . . ,Zl), 



(3.6) 
(3.7) 



(3.8) 

where 
(3.9) 



a(u) — 



hx{u-2r ] ) 
On («) ' 



(3a(u) 



0ii (u - A - 2ar;)0ii(277) 
0ii(u)0ii(A + 2a?7) 
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Let us recall the definition of iV-cycle with step k in Q| . 

Definition 3.1. Let f be a function on and C be a subset of {(x\ + 
miK, . . . , xn + m^jK) | (mi, . . . , tojv) € Z* N } for a certain [x\, . . . , xjv). We call 
C an N -cycle with step k for ^ if 

tec tec 

for any n £ I* N . 

Functions F l (t) and <I>(i) are defined as solutions of the ordinary difference equa- 
tions: 



(3.10) 



(3.11) 



F\t + K) 



a l (t + k) 



F l (t), 



a(i + k) 

We shall give an explicit solution to these equations later in terms of infinite prod- 
ucts. 

The main theorem of this paper is: 

Theorem 3.2. Let 

(3.12) 

M N 

= e°^ u ■ n *i) • n n pin ^ - *o, 

l<i<j<M j=ln=l 

where t = (t\, . . . , tw). Then 

(3.13) f(z)=J2<p(z\t)y(t) 

tec 



is a solution of (2/4), where 
(3.14) 



r-l 



*(f) = e^^Ba+^xfazN, ... , Zl )B a+2 ^ 2 (t 2 ; z N , 



,zi)--- 



a=0 



■ ■ ■ B a+M ^ a _ M (t M ; z N , . . . , zi)Q a _ M (z) 

and C is an N -cycle with step k for <p(z\t )^(t ). 

Graphically, a summand in (3.14) is denoted as in Figure [l4|. 

Remark 3.3. The vector ^(t) gives an eigenvector of the transfer matrix of the 
XYZ typ e spin chains when {v,t\,... , tu} satisfies so-called Bethe equations. (Sec 

m, ||), 03) 

The proof of the theorem is essentially the same as that of Theorem 1.4 of §§. 
We first reduce the equation (2.4) including the operator Aj ((2.1) or (2.2)) to the 
system 

(3.15) j (zi -Zj-K)... Rjll'J (zj-i - zj - K)f{zj) = 
using the unitarity (1.11). 
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a - 1 J a - 2 ■■■ a-M+1 \ a-M 











a-M 


zi -+ 




to 


L M 


a-M 
a-M 


Z 3 * 








a-M 
a + M 


Z N * 








a + M 



+h-l) 



IN 



a 1 a + 1 1 a + 2 ■ ■ ■ a + M - 1 1 a + M 
i i i 



Figure 14. Ba+i,o-i(*i)-"-Ba+M,o-M(*Af)fio-Af(^) 



Assuming the form of the solution ( 3.13 ), we can show that (3.15) is equivalent 

to 

(3.16) 



E 



n /jv,... ,h,lj i, _ 
'^a+JVLa-M^Mj Zn, 



,Zl,Zj 



lj,lx,.. 

,-M 



E< 



^(i1f)^( e 2 *^;:f (t i; . . . , z x ) 



' ' '^a+MM-MV-M'iZ^ZN, 



(Zl, . . . ,2jv,2j-)). 

Zj+i,^--!, . . . ,Z\,Zj). Equiva- 



Here, for example, {zn, ■ ■ ■ ,Zi, Zj) mea ns fa y, . 
lence of the left hand side of (3.15) and ( 3.16 ) is proved as illustrated in Figure [L5| 
We move the line with the spectral parameter Zj + k, repeatedly using the pro- 
cedures Figure |J, and then Figure |[ Because of the admissibility condition (1.5) 
and the equation (1.14), the IRF-type weights appearing in Figure |l5| (crossings of 
dashed lines) are equal to the unity. The right hand side of (3.16) is derived from 
the right hand side of (3.15) in the same manner. 



DIFFERENCE EQUATIONS AND BETHE VECTORS 



15 



a ■ • ■ \ a — M 













Z\ 








Zj-! 






Zj + K 






Z N 







a | • • ■ J a + M 



a I ■ ■ ■ \ a — M 











Zj + K 


t X ••• 




^ ~l 
















1 


ZN 









a 1 ••• 1 a + M 

i i 



Figure 15. Equivalence of (3.15) and ( |3.16|) . 



In the next step we need a counterpart of Lemma 2.3 of 
Lemma 3.4. ( Two-side formula) 
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Fix n (1 ^ n ^ N — 1) and m (1 ^ m ^ M) and Zet a 6e aw integer, ao = a and 
dj =a + 2(h H h ij)r7 /or all j E {1,... , N}. Then 

(3.17) 

AT n 

e n «(*<»**') n n ^-^nn^^-^ 

{1,... ,m}=IUIIieI,i'eII i£l fe=n+l i'ellfc=l 

^a" -J(I)+tl(II) + l,ao+m-l(^i > Z ™> ' ' ' ' Z l) " ' 

' ' ' ^o"'+S(n),oo+lt(n)(**j(i) ' z ™> ■ ■ ■ ' ^^ao'+DCn)^ 1 ' ' • ■ > z «) ® 

' ' ' B oX-it(I)^n-|»(I) ^'«(II) ' 2:JV ' ■ ' ' ' Z "+ 1 ^al + -»(lV i ™ ■ ■ ' ' Z ^)' 

where {1, ... , to} = I U II is a partition, jj denotes the number of elements, i m 
designates an element of I and an element of II. 

Due to (3.6), f/ie n<?/ii /land side of ( 3.17 ) is well-defined. 

One can prove this lemma by induction on m. Using (3.7), (3.8) and a formula 
B l a ^r' h (u;z N ,... ,zi) = 

= A a,b" ' " + 1 ( U '> Z N,--- , z n+l)B l b ^y ' 1 (u; Z n ,... , Z X ) + 

+ B l ^-- ln+1 (u;z N ,... ,z n+ i)D l b ™ a y ll (u;z n ,... ,zi) 

for any integer 6, which is easily derived from (3.2), one can apply the formula (5.6) 
of [IT] to thi s cas e. 

Applying fl3.17p for (m,n) = (M, 1) and (M,N - 1) to the equations fl3.16| ), 
using (1.28) and comparing the coefficients of 

2 m (a+2UU)-2l ] )^ B h r u +K . Z . + K \... 

e - D a+l,a+2tl(II)-2( i -lV^i ' 3 ^1 

' ' ' jS a+t)(II),a+tl(II)-2^ (*»Jcii) + K ' Zj ' + K )^a+|t(II) ( Z 3 + K ) ® 

59 ' CS a+2tl(II)-2i3+l,o-l V il! ^ * " ' ' Z J+1> *J-1) • • • ' Z U'" 

" ' - D a+Af-2i i +)J(II),a-)J(I)V'' i tt(i) ' ZjV ■ • ■ > *j+l > %'-l> • • ■ 

Si a+Af-2i J +)J(II) V Zl ■ • ■ ' Z 3+li ■ • • I *ATj 

of the resulting equation for a partition {1, . . . , M} = I U II, we can show that 
( |3.13j ) gives a solution to (2.4) if tp(z\t) satisfies 

(3.18) 

N 

(p(zj \u + Ks ilU ) n a ^ + K . t ')ii n ah & + k ~ Xfc ) n §h ^ ~ ^ - «) = 



iei.i'eii t'eiifc=i iei 



iGl.i'ell iel i'6llfc=l 

Here in the left hand side 5^ n = if i S I and = 1 if i G II. The function defined 
by (3.12) is a solution of (3.18). Thus we have proved the theorem. 
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Ordinary difference equations and cycles. We now return to the equations 
(3.10) and (3.11) and give the solution to them in the form of infinite products. 

First let us consider equation (3.10). Using the infinite product formula of the 
theta functions (see, e.g., |Mu|), we have 

(3.19) 

6 l (t) =e _ 4mlri (e-^+^jpJoo {pe^-^'-p)^ 



a l (t + k) (q- 1 e- 2wit -' i7Ti ^;p) 00 (pqe 2 ™ t+i7rilr i;p) oc ' 

where p — exp(27rir), q = exp(27riK) and (x;p)oo — [|^_ (1 — p rn x). It is easy to 
see from (3.19) that 

(3.20) 

iqe -2.u + imi v . pj q)oQ ^mt+tiriir,.^ q)oQ 



F l (t) — e ~ i '* il 'n t /K 



satisfies (3.10), at least formally, where (x;p, q)^ = Ilm=oIIn°=o(l —p m Q n x) 
Similarly, equation (3.11) has a formal solution 

(3.21) 

$(t) = e ^ (e -arit ;p)a)(pe a»ft ;j)) 



( qe -2irit-ini V . pj g)oo ( pqe 2^t-4^. p ^ ^ 



(e- 2wit + 4 ^;p,q) 00 (pe 2mt + 4 ^;p,q) oc ' 

In fact, the infinite products in ( |3.20 ) and (3.21) gives meromorphic functions 
on the whole complex plane, provided that Iuik > 0. The function F l (t) has zeros 
at 

(3.22) {t = n + ui\T + m%K + 2lrj \ n G Z, mi G Z>o, TO2 G Z>o} U 

U {t = n — myr — rri2K — 2lr\ \ n G Z, mi G Z>o, mi G Z>o}, 

and poles at 

(3.23) {£ = n + miT + m,2K — 21-q \ n G Z, mi G Z>o, j«2 G Z>o} U 

U {£ = n — m\T — 1712K + 2lr/ \ n G Z, mi G Z>o, m2 G Z>o}. 

Taking these properties of F l (t) into account, we can choose an iV-cycle in ( |3.13 ) 
so that the sum over t G C reduces to a finite sum, if k satisfies a certain ratio- 
nality condition. Suppose that there exist integers (n, mo, mi) G Z> x Z x Z> 
satisfying 

(3.24) nn = mo + miT + AljTj 

for a certain j G {1, . . . , N}. Then all but finite points in the set 

(3.25) C := {t = m° + m?r + inn + 21^ \mEZ} 

fall into the set of zeros of F l (t) (3.22) for any integers m[J and m\ ^ 0. We can 
choose a cycle C so that t n — xj G C for all n — 1, . . . , M and, thanks to the zeros 
of functio n F l j , the sum over t G C in j3.13] ) is essentially a finite sum. Hence 
equation (|3.13| ) gives an analytic solution of the difference equation (2.4). 

4. Comments 



In recent years several difference equations with elliptic coefficients related to 
the g-Knizhnik-Zamolodchikov equations have been proposed: 
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Etingoff's equation: Etingoff @ showed that a (modified) trace of certain in- 
tertwining operators of representations of quantum affine universal enveloping 

algebras, U q (g) satisfies a difference equation with elliptic coefficients. 

Jimbo-Miwa-Nakayashiki's equation: Jimbo, Miwa and Nakayashiki | JMN| 
found a difference equation which should be satisfied by correlation functions 
of eight vertex models. 



Unfortunately we do not know the relation of our difference equation (2.4) with 
any one of above equations. It might be possible that one system turns into another 
by specializ ation of parameters. We can also expect that a quasi-classical limit of 
our system (2.4) is related to the Knizhnik-Zamolodchikov-Bernard equation. 

It is a challenging problem to give a representation theoretical interpretation to 
the so lution ( 3.13 ) like that of the integral solution of the (q-)KZ equations [Ma|, 

ffrJ. 



Appendix A. Review of the Sklyanin algebra 

In this appendix we recall several facts on the Skly anin algebra and its represen- 
tations from [SI | and We use notations in [Mu| for theta functions: 



(A.I) 



ti£Z 



O a b(z;r) = J^exp ( iri +nj r + 2iri(^ 



a) I - + z 



where r is a complex number such that Imr > 0. The Pauli matrices are defined 
as usual: 

<A ' 2) °° - (i !) • - (! ;) ■ °*= (° o') . «■ - G -i 

The Sklyanin algebra, U Ttfl (sl(2)) is generated by four generators S°, S , S 2 , S 3 , 
satisfying the following relations: 

(A. 3) Rniu - v)L 13 (u)L 23 (v) = L 2 3{v)Li 3 (u)R 12 (u - v). 

Here u, v are complex parameters, the L operator, L(u), is defined by 



(A.4) 
where 



a=0 



W„ L (u) 



L(u)=Y,Wt(u)o- a ®S\ 
?u(u;t) 



26 ll (2i r ,r)9 11 (v,r) 



26 11 (27 r ,T)9 W (7 T ,T) 

L , ., 9 01 (u;t) 



26 11 {2r l -T)e m {mT) 1 " a ^ J 2e 11 {2r 1 -T)6 Q1 (r l -T) 
R{u) = R(u; t) is Baxter's R matrix defined by 
(A.5) 

3 

R(u)=Y,W*{u)o- a ®o- a , W*(u) ■.^e ll {2r l ;T)W^{u + i 1 ). 



a=0 
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and indices {0,1,2} denote the spaces on which operators act non-trivially: for 
example, 

3 3 

R 12 (u) = ^W*{u)l®a a ®(T a , L 13 (u) = J2 W a(u)a a ®l®S a . 



The above relation (A.3) contains u and v as parameters, but the commutation 
relations among S a (a = 0, . . . , 3) do not depend on them: 

(A.6) [S a , S }- = -iJaAS 13 , S 7 ]+, [S a , = i[S°, 

where (a,f3,j) stands for any cyclic permutation of (1,2,3), [A, B]± = AB ± BA, 
and J a jj = - Wp)/{W^ - Wq) depend on r and n but not on u. 

The spin I representation of the Sklyanin algebra, p l : U T ^ v (sl(2)) — > Endc(6(jo) 
is defined as follows: The representation space V 1 is 

(A.7) 

V 1 = 6$+ := {f(y) | f(y + 1) = f(-y) = f(y), f(y + r) = exp- 4 ^ 2 ^> /(„)}. 

It is easy to see that diml/' = 21 + 1. The generators of the algebra act on this 
space as difference operators: 

(a.8) ( P l (s°)f)(y) = s ^y- l ^y + *)) «.(-» l ^)f(y »?) 



^i(2y;r) 
where 

So(y) = 11 (r];T)d 11 (2y,T), si(y) = d 10 {n;T)9 10 (2y;T), 
s 2 (y) ^i8 Q0 {'n;T)e QQ (2y;T), s 3 (y) = d i(r); t)6 01 (2y; t). 

These representations reduce to the usual spin I representations of U(sl(2)) for 
Ja(i —* (n — > 0). In particular, in the case I = 1/2, S a are expressed by the Pauli 
matrices a a : Let us identify 6qq and C 2 by 

9 Q0 (2y;2r)~9 10 (2y;2T)^ (J 

(A - 9) )o\ 
M2y; 2r) + M2y; 2r) ^ (J J 

Under this identification S a have matrix forms 
(A.10) p 1 / 2 (S a )=9 11 (2r,;T)a a . 
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